Complex dislocation dynamics in ice: experiments 
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We present a statistical analysis of the acoustic emissions induced by dislocation motion 
during the creep of ice single crystals. The recorded acoustic waves provide an indirect 
measure of the inelastic energy dissipated during dislocation motion. Compression and 
torsion creep experiments indicate that viscoplastic deformation, even in the steady-state 
(secondary creep) , is a complex and inhomogeneous process characterized by avalanches in 
the motion of dislocations. The distribution of avalanche sizes, identified with the acoustic 
wave amplitude (or the acoustic wave energy), is found to follow a power law with a cutoff 
at large amplitudes which depends on the creep stage (primary, secondary, tertiary). These 
results suggest that viscoplastic deformation in ice and possibly in other materials could 
be described in the framework of non-equilibrium critical phenomena. 
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I. INTRODUCTION 

ithin a material-dependent range of temperature 
and stress |^], the (visco)plastic deformation of 
crystalline materials involves the motion of a large 
number of dislocations. The importance of collec- 
tive effects, through clastic interactions, on disloca- 
tion dynamics and the heterogeneous character of 
plasticity has been recognized for a long time (see 
e.g. 1^). Only relatively recently, 2D ||,^ as well 
as 3D [g| numerical simulations of collective disloca- 
tion interactions have revealed a spontaneous self- 
organization of dislocation patterns with walls and 
cells, in agreement with observations |^ . How- 
ever, collective dislocation dynamics has received 
much less attention, with the exception of the spe- 
cial case of Portevin-Le Chatelier (PLC) effect (see 
below). This may be due to the difficulty to moni- 
tor experimentally the variations of dislocation mo- 
tion in time, energy and space domains. Here, we 
report an experimental analysis of acoustic emis- 
sions (AE) generated by dislocation motions which 
allow to reveal the intermittent and jerky charac- 
ter of collective dislocation dynamics. Experiments 
are performed on almost disorder free (except dis- 
locations) ice single crystals. In a companion paper 
|7j, we present numerical simulations which repro- 
duce the main statistical characteristics of the AE 
measurements. This work (experiments and model) 
show that complex and jerky plastic flow is not re- 
stricted to the PLC effect, and rather the rule than 
the exception. 



II. ACOUSTIC EMISSION TO MONITOR 
DISLOCATION DYNAMICS 

In solid materials, sudden local changes of inelas- 
tic strain generate AE waves ||] . The sources can be 
crack nucleation and propagation, twinning, or dis- 
location motion. In the present work, AE is used to 
monitor dislocation dynamics. In our experiments, 
given the amplitude threshold and the frequency 
range of our transducer, the detected AE are un- 
likely to be the result of a single moving disloca- 
tion but most probably are related to synchronized 
accelerations of dislocations, called plastic instabil- 
ities or dislocation avalanches From the theo- 
retical analysis of Rouby et al. | Iffl one can relate 
the maximum amplitude of the acoustic wave A re- 
sulting from a plastic instability, to the number of 
involved dislocations n and their velocity v p|: 



A = k 



nLhvtf) 



(1) 



where fc is a coefficient related to material prop- 
erties and the piezoelectric constant of the trans- 
ducer (dimension: V/m), b is the Burger s vec- 
tor, L is the length of the n moving dislocations, 
to is the travel time of the acoustic wave through 
the transducer (considered to be constant), and d 
is the source/transducer distance (supposed to be 
large compared to L). In this crude model, L and 
V are supposed to be identical for all the involved 
moving dislocations. The dislocation velocity v is 
considered to be zero before and after the event and 
constant during the event. The term \/d represents 



the geometrical attenuation of the acoustic waves. 
From (1), Weiss and Grasso Q showed that the AE 
ampHtude ^ is a measure of the local strain associ- 
ated with the dislocation avalanche. Therefore, the 
rate of global AE activity is an indirect measure of 
the global strain rate of the sample. 

The AE energy radiated by the acoustic wave is 
proportional to A^. According to Kiesewitter and 
Schiller jll|, the energy dissipated by viscoplastic 
deformation during an event also scales with A^. 
This results from an expression given by Eshelby 
|l2| for the energy dissipated at the source by a 
single screw dislocation of length L moving at a 
velocity v: 



E = KL^V 



(2) 



where K \s a. coefficient depending on material con- 
stants, including the shear modulus and the veloc- 
ity of acoustic transverse waves, and h the Burger's 
vector. A comparison of (1) and (2) with n = 1 
shows that E ^ A^. This scaling has been observed 
during our experiments. Therefore AE allows to 
study during deformation the dislocation dynamics 
in energy, time, and possibly space (if AE sources 
locations are determined with the help of multiple 
transducers) domains. In the present study, only 
two transducers were used, which did not allow 3D 
localization. 



III. EXPERIMENTAL PROCEDURE 

As a model material to study dislocation dynam- 
ics from acoustic emission, ice provides the follow- 
ing advantages: single crystals or polycrystals 
with various microstructures can be easily grown 
in the laboratory; transparency allows verification 
that AE activity is not related to microcracking; 
and an excellent coupling between the ice and the 
AE transducers can be obtained by fusion/freezing. 
Within the range of temperature and stress cor- 
responding to our experimental conditions, diffu- 
sional flow is not a significant mechanism of de- 
formation in ice, and viscoplastic deformation oc- 
curs by dislocation motion Hexagonal ice Ih 
presents a very strong plastic anisotropy of the sin- 
gle crystal ]14|: viscoplastic deformation of single 
crystals occurs essentially by basal glide. 

Uniaxial compression creep as well as torsion 
creep experiments were performed at -10 oC on ar- 
tificial ice single crystals (150 mm X 75mm), each 
of them constituted by several steps of constant ap- 
plied stress. These experiments were similar to ex- 
periments performed previously and which are de- 
tailed in and For the present new set of 
experiments, a new AE recording device was used 
which allowed to record, for each event, different 
characteristics including the arrival time (at a pre- 
cision of 0.1 fis), the amplitude, the acoustic energy 



or the average frequency. Owing to a greater sensi- 
tivity of the recording device, much better statistics 
on AE events were obtained compared to previous 
work . The frequency band width of the trans- 
ducers was 10-100 kHz. During an experiment the 
event amplitude threshold was adjusted to 30 dB, 
or 310"^ Volts, i.e. about 5 dB above the noise 
level. The dynamic range between the amplitude 
threshold (310 —3 V) and the maximum recordable 
amplitude (10 V, or 100 dB) was 70 dB, i.e., 3.5 
orders of magnitude. The corresponding dynamic 
range for energies was 7 to 8 orders of magnitude. 

Several independent evidences for recorded AE 
to be related to dislocation motion, including an 
absence of microcracks or a proportionality between 
the global strain-rate and the global AE rate have 
been detailed elsewhere IqIJTHI . 



IV. STATISTICAL ANALYSIS OF AE 
AMPLITUDES AND ENERGIES 

From equation (1), the acoustic activity A{t) 
monitor the evolution of dislocation activity and 
viscoplastic dissipation. Figure 1 shows this evo- 
lution during the first loading step of a compres- 
sive creep test. It appears immediately that the 
instantaneaous viscoplastic dissipation is intermit- 
tent and jerky, showing bursts of activity that can 
be considered as avalanches of dislocations moving 
collectively in the material. Although the cumu- 
lated acoustic activity, i.e. an estimate of the strain, 
appears relatively smooth at the time scale of sev- 
eral thousands of seconds (Figure 1), At smaller 
time scales, or in terms of dissipated energy E{t), 
this avalanche behavior is more obvious (Figure 2): 
dislocation dynamics during viscoplastic deforma- 
tion is heterogeneous in time and energy domains. 
This recalls the jerky flow associated to Portevin- 
Le Chatelier effect occuring in different alloys (see 
e.g. ||l^,|l^). However, whereas PLC effect is at- 
tributed to the dynamic interaction of two defect 
populations, namely mobile dislocations and solute 
atoms [^1, the avalanche behavior reported here 
results on y from collective dislocation interactions 
(solid solutions in the artificial single crystals are 
below 10-10 mole frac. citelS ). For this reason, 
we believe that this kind of behavior should not be 
restricted to ice, as confirmed by a numerical mod- 
elling of this effect in a companion paper . 

This avalanche behavior is recovered whatever 
the type of loading (compression or torsion) , the ap- 
plied stress, or the creep stage (primary, secondary, 
tertiary). However, interestingly and unexpectedly, 
the classical distinction between these creep stages 
corresponds also to different dislocation dynamics: 
primary and tertiary c reep are associated to a pro- 
nounced jerky character and some very large events, 
whereas secondary creep appears smoother, at least 



2 



at a large time scale. This difference, visible on fig- 
ure 1, is clearer on the amplitude distributions. Sec- 
ondary creep is characterized by power law statis- 
tics with an exponential cutoff at large amplitudes 
(or energies) (Figure 3): 

N{A) - yl-^e-^/^'= (3) 

or the discrete distribution calculated with linear 
bins, or 

N{A' > A)^ A-^+^e-"^/^" (4) 

for the cumulative distribution. Because of the scal- 
ing E A^, similar relations are obtained for en- 
ergy distributions with an exponent (r -I- l)/2. For 
secondary creep under compression, the exponent r 
of the discrete distributions, estimated from a lin- 
ear fit on a log-log plot below the amplitude cutoff 
Ac, is found to be stable through creep deforma- 
tion, independent of the applied stress and equal to 
2.0 ± 0.05. Similarly, Ac is found to be constant 
during a given loading step (Figure 3) and does not 
seem to vary significantly with the applied stress 
around Ac ~ 0.4 — 0.5V. This means that power 
law scaling is clearly observed over 1.5 orders of 
magnitude, i.e. 3 to 4 orders for energies. The 
stability of this scaling justifies, on the dynamical 
point of view, the term steady-state sometimes used 
to identify secondary creep. 

On the other hand, statistics of dislocation 
avalanches change with the stage of creep. Pri- 
mary as well as tertiary creep are characterized by 
a power law scaling without detectable upper cut- 
off over the experimental scale range {A ~ lOT^). 
In other words, A^ if any, is rejected towards much 
larger scales compared to secondary creep. This 
transition from secondary to tertiary creep is visi- 
ble on figure 4. Slightly larger r exponents, around 
2.15 ± 0.05 seem also to be associated with primary 
and tertiary creep. The variation of AC with the 
creep stage suggests that this upper cutoff is not 
a finite-size effect. Under torsion creep, power law 
scaling was also observed without detectable upper 
cutoff, but with a smaller exponent, r = 1.7. This 
difference of t on the loading mode is not explained 
so far. 



V. TIME PATTERNING 

To study the time patterning of dislocation dy- 
namics, we calculated the correlation integral C{t): 

C{t) = J^^fi^t < t) (5) 

where N is the total number of events and n{/S.t < t) 
is the number of pairs of events (not necesseraly suc- 
cessive events) separated by less than a time t. C{t) 



is therefore simply the probability for two acoustic 
events to be separated in time by less than t. Note 
that this analysis does not take into account the am- 
plitude of the events. For a random Poisson process 
with uncorrelated (identical) events, C{t) ~ A 
typical profile of C{t) is given on figure 5. Whereas 
a scaling C{t) ~ t is observed above a threshold 
tc, the events are clustered below tc (the proba- 
bility C{t) is larger than expected for uncorrelated 
events). This behavior was observed whatever the 
type of loading (compression/torsion), the apphed 
stress or the creep stage, but tc was found to follow 
the evolution of < At >, the average time interval 
between two successive events. These observations, 
along with detailed examinations of the data files, 
can be explained as follows: dislocation avalanches 
consist of a mainshock correlated in time with a se- 
quence of few aftershocks. Note that these shocks 
are distinct dynamical events characterized by dis- 
tinct acoustic waves. At larger time scales, vis- 
coplastic deformation consists of avalanches which 
are themselves uncorrelated in time. When the AE 
activity is very high, the frequency of avalanches is 
so high such that the aftershock sequences are not 
revealed by the correlation integral analysis. 

VI. DISCUSSION AND CONCLUSION 

From a statistical analysis of acoustic emis- 
sions generated by dislocation motion, we have 
shown that: (i) viscoplastic deformation during 
the creep of ice occurs by a succession of dislo- 
cation avalanches. This avalanche behavior, in 
very pure ice single crystals, results only from 
dislocation-dislocation elastic interactions, (ii) dis- 
location avalanches consist of a mainshock corre- 
lated in time with a sequence of few aftershocks. 
These avalanches are themselves uncorrelated in 
time, (iii) the distributions of AE amplitudes (the 
"size" of the avalanche) or AE energies follow a 
power law scaling with an exponential upper cut- 
off which depends on the creep stage. For a specific 
creep stage, the statistical properties are stable and 
do not vary with time or the applied stress. Such 
power law scaling of dislocation avalanches suggests 
that viscoplastic deformation in ice and possibly in 
other materials could be described in the frame- 
work of non-equilibrium critical phenomena . In 
particular, driven-dissipative systems have the ten- 
dency to develop singular response functions and 
avalanche like behavior depending on the driving 
mechanism. This critical behavior is, in fact, re- 
sulting from the collective behavior of the many 
microscopic degrees of freedom of the system re- 
sponding to the external perturbation. It is then 
of a major interest to develop theoretical models 
that can relate the observed experimental behav- 
ior with the microscopic dynamics of dislocations 
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and its interactions with the external driving fields. 
In the companion paper [7] we provide an attempt 
in this direction by developing a numerical model 
that successfully reproduces the avalanche behav- 
ior of dislocation motion, (iv) whatever the very 
nature of the involved critical phenomena, Weiss 
et al. have shown that power law distributions of 
dislocation avalanches imply that large plastic in- 
stabilities account for most of the viscoplastic de- 
formation, rather than independent movements of 
individual dislocations. 

The results presented in a companion paper [7] 
represent a further evidence to argue that this 
avalanche behavior could not be specific to ice, 
rather a common feature of collective dislocation 
dynamics. 
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FIG. 1. Instantaneous as well as cumulated acous- 
tic activity (amplitude) during the first loading step 
of a compression creep test (applied stress=0.58 MPa). 
Transition from primary to secondary creep. 
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FIG. 2. Detail of figure 1 within a time window: Cu- 
mulated amplitude (solid line) as well as cumulated en- 
ergy (dashed line) 
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FIG. 3. Cumulative 

amplitude distributions N(A' > A), normalized by the 
total number of events in the time window, A^, during 
secondary creep under compression (fourth step of load- 
ing; applied stress=1.29 MPa). The three sets of data 
(2500 events each, or a duration of about 110-120s) cor- 
respond to three non-overlapping time windows. The 
dashed line corresponds to r = 2 (r — 1 = 1 for the 
cumulative distribution) 



FIG. 4. Cumulative amplitude distributions 
N{A' > A), normalized by the total number of events 
in the time window. A'', during the sixth step of loading 
of a compression creep test (applied stress=1.64 MPa): 
Transition from secondary to tertiary creep. The three 
sets of data correspond to three non-overlapping time 
windows of duration 258s. This duration corresponds to 
about 5000 events for two windows ( and )) and about 
71000 events for the last one (F, tertiary creep). 
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FIG. 5. Correlation integral C{t) calculated on two 
non-overlapping time windows (2500 events each) dur- 
ing the second loading step of a compression creep test 
(applied stress=0.71 MPa). The dashed line corre- 
sponds to C{t) ^ t. The sharp decrease of C{t) below 
510^^s results from the duration of the longuest events: 
the experimental device is unable to detect a new event 
as long as the former event is not finished. This limits 
artificially C{t). 



